Temperature dependent change in the symmetry of the order parameter in an 
electron-doped high-temperature superconductor 



O 

o 

(N 



o 
o 

I 

Oh 

^ , 



I 

O 

o 



> 

(N 
(N 
O 

o 



I 

o 
o 



Hamza Balci and R. L. Greene 
Center for Superconductivity Research, Department of Physics, University of Maryland, College Park, MD-20742 

(February 2, 2008) 

We present specific fieat measurements wliicli sliow an unexpected pfiase transition from d-wave 
symmetry to s-wave symmetry as tlie temperature is reduced in electron-doped Pr2-a;Cei:Cu04_i 
(PCCO), in both optimal and over-doped single crystals. The field dependence of electronic specific 
heat (Ce;) is linear at T=2K, consistent with s-wave symmetry, and non-linear, consistent with 
d-wave symmetry, at T>3K. This behavior is most consistent with a phase transition in the sym- 
metry of the order parameter as the temperature is changed. Such a phase transition could be an 
explanation for the previous controversial results of different experiments, performed at different 
temperatures in the electron-doped cuprates. 



Particle-hole symmetry is an important ingredient of 
any theory attempting to explain the mechanism of high 
temperature superconductivity. Despite its fundamen- 
tal importance, this issue has not been resolved yet, i.e. 
there is no consensus as to whether the electron-doped 
(n-type) and hole-doped (p-type) cuprates are essentially 
the same or very different. The symmetry of the order 
parameter is one of the important parameters to com- 
pare the n-type cuprates to the p-type family. Unlike 
the p-type cuprates, which are now generally considered 
to have d-wave symmetry [1,2], the pairing (gap) sym- 
metry of the n-type cuprates continues to remain contro- 
versial, despite the enormous improvement in the sample 
quality and experimental techniques since their discovery 
over 10 years ago. In addition to experiments that sup- 
port an s-wave symmetry [3] or a d-wave symmetry [4-8] , 
there are a few experiments that suggest a transition from 
d-wave symmetry to s-wave symmetry as the doping is 
increased [9,10]. 

Specific heat has traditionally been used to study the 
superconducting (SC) gap since it is a direct measure- 
ment of the density of states (DOS) . The difficulty of de- 
tecting the small T^-type temperature dependence ( 5% 
of the phonon specific heat at 4 K in YBCO [12]) ex- 
pected at zero field from a d-wave superconductor [11] 
makes field dependence studies the preferred method of 
studying the gap symmetry. In a simple isotropic s-wave 
superconductor, the specific heat is linear in magnetic 
field for fields larger than Hd (II>Hci) [13]. For a d-wave 
superconductor, Ce; has a y/H-type field dependence in 
the low temperature (T << Tc) and intermediate field 
regime (i?ci « H « Hc2) [14]. 

In our case measuring the field dependence of Gei is 
particularly advantageous since our PCCO crystals do 
not have electronic or nuclear Schottky contributions (see 
Fig. 1) in the temperature (T>2K) and field (H<10T) 
range of our measurements. The field independent ad- 
denda leaves Cei the only field dependent part of in our 
specific heat measurements. Therefore, our field depen- 
dent data was acquired at a constant temperature while 
the magnetic field is ramped up (zero field cooled) or 



down. No hysteresis has been observed in our specific 
heat between increasing and decreasing field. 

We measured the heat capacity of Pr2_a;Cea;Cu04_5 
(PCCO) single crystals grown by the directional solidifi- 
cation technique with cerium dopings of x=0.150±0.005 
and 0.165±0.005. The non-superconducitng as-grown 
crystals are annealed at 900 °C (in an inert atmosphere 
of flowing argon gas) in order to attain superconductiv- 
ity. The cerium concentration of the crystals was deter- 
mined using WDX to an accuracy of ±0.005. The ex- 
periments were performed in a Quantum Design PPMS 
with a modified sample holder using thermal relaxation 
calorimetry [15]. Au-7%Cu wires (1-3 mil diameter) are 
used as a weak link to connect the holder to the ther- 
mal bath. The thermometers were calibrated at different 
magnetic fields in order to take into account their small 
magnetoresistance. Our setup was tested by measuring a 
3 mg high purity copper sample, and a 3.2 mg Nb sample, 
and our measurements were within 5 % of the standard 
values. 

Fig. 1 shows temperature dependence data on an 
optimally-doped PCCO crystal with Tc=22±2 K and 
mass of 3.2 mg in several magnetic fields (II//c-axis). 
The absence of any Schottky upturn at low temperatures 
shows that the sample is free of any measurable magnetic 
impurities. The specific heat in this case can be repre- 
sented as C = [7(0) + 7(i?)]T + /3T3, where 7(0) is the 
coefficient of the zero field residual heat capacity, ^{H) 
is the coefficient of the field dependent Ce/, and (3 is the 
coefficient of the phonon specific heat. The gap sym- 
metry is investigated by studying j{H). For H > Hc2, 
7(_ff) — 7„, where 7„ is the Sommerfeld constant. The 
7(0) term is also found in hole-doped cuprates, and there 
are several proposals to explain its origin [16-18]. By 
making a linear fit to the C/T vs T^ data shown in the 
inset of Fig. 1 {(3 is kept constant for the two fields), we 
find 7(0)=2.1±0.2 mJ/moleK2, 7„=5.3±0.2 mJ/moleK^, 
and /3=0.23±0.02 mJ/moleK^, which results in a Debye 
temperature 6'd=390±15 K. 

The Debye temperature is the same (within the accu- 
racy of our measurements) in all the crystals we studied. 
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and the value of 7(0) changes between 1.4-2.2 mJ/moleK^ 
in the different optimaUy-doped crystals we studied. The 
7„ for different optimally-doped crystals depends on the 
annealing time (oxygen content) and we have observed a 
variation between 4.0-5.3 mJ/moleK^. The upper criti- 
cal field Hc2 determined from specific heat measurements 
also depends on the annealing time. Fig. 1 shows that 
the upper critical field for this sample is B.c2 <8T since 
there is no difference in the specific heat data between 
8T and lOT. The Bc2 of the different optimally-doped 
samples changed between 5.0-7.0 (±0.5T) depending on 
the annealing time of the crystals. 

Fig. 2 shows the field dependence of an optimally- 
doped crystal (Tc=23±3K, mass=3.4 mg). Only the 
field dependent part of Cgi is shown in Fig. 2 since the 
zero- field specific heat, which includes the phonon con- 
tribution (/JT"*) and the zero field residual heat capac- 
ity (7(0)T), is subtracted out. The field dependence is 
clearly non-linear at T>3K. However, a dramatic change 
in the field dependence is observed when the temperature 
is reduced to 2K, where a linear field dependence dom- 
inates up to about 3T. Fig. 2-B shows that the change 
in the field dependence of Cei between 4.5K and 2K is 
accompanied by a suppression of the electronic 
specific heat. This suggests that the density of elec- 
tronic states is suppressed, or possibly a gap is opened 
at the Fermi level, since the electronic specific heat is 
a measure of DOS. 

This change in the field dependence of Get from lin- 
ear at T=2K to non-linear at higher temperatures has 
been observed on four optimally- cerium doped crystals 
with different oxygen concentration. In addition to 
the optimally-cerium doped crystals, crystals of higher 
cerium doping also show a similar field dependence. 
Fig. 3-A shows the field dependence of specific heat 
for an over-doped crystal (x=0.165±0.005, Tc=15±1.5K, 
mass=2.2 mg). The trend from linear to non-linear field 
dependence as the temperature is increased is also ev- 
ident in this crystal. In addition, Fig. 3-B shows that 
there is a small suppression of Gei, and hence electronic 
DOS, as the temperature is reduced from 3K to 2K. 

There arc several possible scenarios to explain such an 
unusual change in the field dependence of Cei- a phase 
transition in the symmetry of the order parameter from 
d-wave at high temperatures (T>3K) to s-wave at lower 
temperatures(T<2K); an anisotropic s-wave gap which 
has a minimum around 3K; or a vortex- vortex interaction 
effect which results in a non-linear field dependence in 
an s-wave superconductor at T>3K. Now we will discuss 
each of these scenarios in more detail. 

Vortices can interact with each other via quasiparti- 
cle (QP) transfer between their cores. At high mag- 
netic fields (large number of vortices) or high temper- 
atures (larger vortex-core size) the quasiparticle wave- 
functions in the core of one vortex overlap with the 
quasiparticle wavefunctions in neighboring vortices, and 
hence inter-vortex quasiparticle transfer becomes possi- 
ble. These inter-vortex QP transfers result in a shrink- 



ing of the vortex cores [19]. The calculations of Ichioka 
et al. [20], which took into account these vortex-lattice 
effects, showed that Cei ocHO *^^ at T=0.5Tc. Calcu- 
lations of Miranovic et al. [21] at lower temperatures, 
where less overlap between vortex-cores reduces these ef- 
fects and the size of the vortex-core is approximately con- 
stant for H<0.5Hc2, showed that Cei c)cH for H<0.5Hc2 
at T=0.1Tc. 

These ideas were experimentally supported by the 
non-linear field dependence of G^i observed in some s- 
wave superconductors, e.g. VsSi [22], NbSe2 [23], and 
CeRu2 [24]. Despite some complications that were dis- 
covered after the original work, a vortex-lattice transfor- 
mation in VaSi [25] and anisotropic gap in NbSe2 [26], 
the fundamental mechanism of vortex-core shrinking due 
to vortex-vortex interaction is still consistent with the 
non-linear field dependence observed in these materi- 
als [25]. 

The change in the field dependence of Cgi in PCCO 
is in qualitative agreement with these theories for an s- 
wave superconductor. However, in PCCO the change 
in the field dependence occurs over a very narrow tem- 
perature range Ce; cxH*^-^ at Tw0.2Tc and Cei ocH at 
TwO.lTc) compared to the gradual and slower change 
based on vortex-lattice effects predicted in Rcf. [20,21] 
for a conventional s-wave superconductor {Cei ocH°'^^ at 
T=0.5Tc and Cei ocH at T=0.1Tc). Hence, it is very 
unlikely that the non-linear field dependence we observe 
in PCCO is due to a vortex lattice effect. Regardless of 
the source of the non-linear field dependence at higher 
temperatures, the linear field dependence that we ob- 
serve at T=2K, which is very similar to what is predicted 
by Ref. [21] for s-wave superconductors, shows that the 
symmetry of the order pareimeter is s-wave at this 
temperature. 

Another possible explanation for our data is an 
anisotropic s-wave gap. However, it is difficult to differen- 
tiate between an anisotropic s-wave superconductor and 
a d-wave superconductor using thermodynamic proper- 
ties because what is usually measured ex perimentally is 
QP excitation energy Ek = \Ji1 + jA^I^, where is the 
single particle energy relative to the Fermi level and 
is the SC gap, which does not depend on the a phase of 
the SC gap. In particular, measurements on YNi2B2C, 
a highly anisotropic s-wave superconductor, have shown 
Cei ocH°'^ very similar to d-wave superconductors [27]. 

For an anisotropic gap the field dependence could be 
linear or non-linear depending on whether the energy 
of the QP excitations is below or above the SC gap 
minimum (Amj„). If we consider our data in Fig. 2, 
^min ~ S/cb («0.26 meV), since Cei is linear in H below 
T=3K, and non-linear at higher temperatures. However, 
there is a problem with this picture. In an anisotropic 
s-wave gap the magnetic field should have a similar effect 
on the QP excitation spectrum as the temperature. The 
QP's are excited due to a magnetic field H by an energy 
E// « AovW^ [28], where Aq is the gap maximum 
(Ao «4meV for PCCO [9]). At T=3K in Fig. 2, the field 
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dependence is non- linear even at fields as low as 0.1 T. 
Yet at T=2K the field dependence remains linear up to 
3T. Due to the similar effects of magnetic field and tem- 
perature on the QP excitation spectrum in an anisotropic 
s-wave gap, for T=2K we would expect the field depen- 
dence to become non-linear at much lower fields than 3T 
(H<1T). Since this is not what we observe in our data, 
we believe that the change from linear field dependence 
to non- linear field dependence in Cgj of PCCO is not due 
to an anisotropic s-wave gap. 

Finally we consider a phase transition in the symme- 
try of the order parameter from s-wave (T<2K) to d- 
wave (T>3K). This interpretation is particularly attrac- 
tive since it has the potential of reconciling the results of 
many confiicting experiments. The SQUID, ARPES, and 
Raman [5-7] measurements that suggested d-wave sym- 
metry in electron-doped cupratcs arc almost exclusively 
performed above 4.2K. On the other hand penetration 
depth (0.4K-Ref. [3]) and point contact tunnelling spec- 
troscopy (1.8K-Ref. [9]), are performed at T<2K and 
they are consistent with s-wave symmetry. The only 
inconsistent experiment with this picture is a penetra- 
tion depth experiment which suggested d-wave symmetry 
down to T=0.4K(Rcf. [4]). 

The above qualitative picture of a phase transition is 
also consistent with quantitative analysis both at 4.5K 
and 2K. The Cei of a d-wave superconductor has a y/H- 
type field dependence given by [29]: 

(o \ 1/2 / TT \ 1/2 

This equation is valid at T < and Hd « «; Hc2- 
Fig. 2-B shows a d-wave fit given by Eq.(l) to the 4.5 K 
data and a linear fit to 2K data. In calculating A, the co- 
efficient of the term, 7„=4.2 mJ/mole K^, Hc2=7T, 
and a=0.7 were used. The 7„ and B.c2 are determined 
from our specific heat measurements. Using these param- 
eters a reasonably good fit of our data to 'y{H) = AH^/^ 
is obtained at T=4.5K (Fig. 2-B). An s-wave supercon- 
ductor on the other hand has a linear field dependence: 
Cei = KjnT H / Hc2{T) , whcrc K is a geometrical constant 
of order 1. Using 7„=4.2 m J /mole and Hc2= 7 T at 
T=2 K we find the theoretical Cgi = 1.2kH. A linear fit 
to our data yields Cei = 2. OH. A comparison with data 
results in «;=1.7, which is a reasonable value for k. 

Such a phase transition should also be observed in the 
temperature dependence of the specific heat. We esti- 
mate the magnitude of this feature to be approximately 
5% of the total heat capacity around 3K (the suppression 
in Fig.2-b is approximately 5% of the total heat capac- 
ity), and expect it to be spread over a temperature range 
of 1-2K. We were not able to resolve this transition in 
the temperature dependence of the specific heat because 
of our limited resolution (5%), however more sensitive 
measurements using ac heat capacity arc in progress. 

The data on the over-doped crystal (shown in Fig. 
3) is also consistent with s-wave theory at T=2K. Us- 



ing IIc2= 3.6 T, 7„=2.8 mJ/mole (determined as for 
the optimally-doped crystal), an s-wave estimate of 1.6k 
is obtained for the slope at T=2K. A linear fit to the 
data yields 2.2 mJ/mole K T. This results in k;=1.4 for 
the over-doped crystal, which is again a reasonable value 
for K. The low Hc2 (Hc2=2.90.2T at T=3K) makes a 
comparison with d-wave theory difficult since the con- 
dition Hci<CH<CHc2 is only satisfied in a very narrow 
field range. However, the trend from linear to non-linear 
field dependence as the temperature is increased suggests 
that the x=0.165 crystal behaves similarly to the x=0.150 
crystals. 

The phase transition suggested by our data is com- 
patible with a theoretical model proposed by Khodel et 
al. [30] which suggests a phase transition in the gap sym- 
metry (from nodal to fully-gapped) in n-type cuprates 
based on the location of the Fermi surface "hot-spots". 
The consistent picture formed by this model and our data 
implies that the hot spots and hence the antiferroniag- 
netic spin fluctuations are an important ingredient of su- 
perconductivity in the n-type cuprates. 

In conclusion, the results of our specific heat measure- 
ments can be summarized as: 

• a dramatic change is observed in the magnetic field 
dependence of the electronic specific heat in a very 
narrow temperature range; non-linear at T>3K 
and linear at T=2K. 

• there is a suppression of DOS, consistent with open- 
ing of a gap at the Fermi surface, as the field de- 
pendence changes from non- linear to linear. 

In light of these observations we conclude that the sym- 
metry of the order parameter in the n-type cuprates is not 
a simple d-wave symmetry as in the p-type compounds. 
In particular the s-wave symmetry observed at T=2K is 
a definite evidence for this difference. The change in the 
field dependence of G^i is most consistent with a phase 
transition in the symmetry of the order parameter from 
d-wave to s-wave. Such a phase transition in the symme- 
try of the order parameter can explain most of the prior 
conflicting symmetry experiments, since almost all ex- 
periments that suggest d-wave symmetry are performed 
at T>4K, and those that suggest s-wave symmetry are 
performed at T<2K. 
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FIG. 2. A-Field dependent part of the electronic spe- 
cific heat (C(II)-C(O)) for an optimally-doped PCCO crystal 
(Tc=23±3K) between 2-4.5K. At 2K the field dependence is 
linear consistent with a fully-gapped order parameter. At 
slightly higher temperatures the field dependence changes 
from linear to non-linear. B-The field dependent specific heat 
is divided by temperature to show the suppression in the spe- 
cific heat, and hence DOS. The curved line is a d-wave fit to 
4.5K, and the dashed line is an s-wave fit to 2K data. 
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FIG. 1. Temperature dependence of specific heat for an op- 
timally-doped crystal, Pri.85Ceo.i5Cu04 (Tc=22±2K), at dif- 
ferent magnetic fields (H//c-axis). The saturation of the spe- 
cific heat, i.e. field-independent, at 8T suggests that IIc2 <8T. 
The inset shows C/T vs T^ for OT and lOT magnetic fields 
from which 7(0), "f„, and /3 are extracted. 
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FIG. 3. A- Field dependence of specific heat for an 
over-doped PCCO crystal (Tc=15±1.5K). The field depen- 
dence of the specific heat is similar to the optimally-doped 
crystals, i.e. linear at 2K and non-linear at higher tempera- 
tures. B-The field dependent specific heat is divided by tem- 
perature to show the suppression in the specific heat, and 
hence DOS, similar to the optimally-doped crystal. 
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